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Levy flights can be described using a Fokker-Planck equation, which involves a fractional 
derivative operator in the position co-ordinate. Such an operator has its natural expression 
in the Fourier domain. Starting with this, we show that the solution of the equation can 
be written as a Hamiltonian path integral. Though this has been realized in the literature, 
the method has not found applications as the path integral appears difficult to evaluate. 
We show that a method in which one integrates over the position co-ordinates first, after 
which integration is performed over the momentum co-ordinates, can be used to evaluate 
several path integrals that are of interest. Using this, we evaluate the propagators for (a) 
free particle (b) particle subjected to a linear potential and (c) harmonic potential. In all 
the three cases, we have obtained results for both overdamped and underdamped cases. 



I. INTRODUCTION 

The field of anomalous transport has gained widespread interest due to its potential to explain 
several phenomena which fall outside the realm of simple Brownian motion [l- 17[. The continuous 



time random walk (CTRW) is a model used to understand anomalous transport 13" 21| where jump 
lengths and waiting time distributions are both drawn from a probability distribution function 
'ip{x,t). The jump length distribution is given by 



X{x) = I dt ij{x,t), (1) 
and the waiting time distribution is given by 

w{t) = I dx ij{x,t). (2) 



Levy flights belong to a special class of CTRW 0, for which, ij{x,t) = \{x)w{t). The 

waiting time distribution, w{t), is narrow, and this makes the process Markovian in nature 24l ]. 
The probability distribution for the jump length, X{x), is Levy stable. The easiest way to specify 
it is to use its Fourier transform. The most general representation of the Fourier transform of a 



Levy stable distribution is given by [2J], 
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and < a < 2 is the Levy index, — 1 < /? < 1 is the skewness parameter, ^ E M is the shift param- 
eter, and a G R"*" determines the strength of the noise. We will write the position representation 
for a Levy stable distribution as Lq,^^(x; cr). This distribution has the property that it can be 
rewritten as 



24| 
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The shift and strength parameters {f^,cr) have been absorbed into just one term in the argument. 
We will adopt an easy to write notation for the Levy stable distribution, which is, ^i^a,/3 (^^^^)- If 
we consider the such a jump length distribution, then the Fourier transform of A(x) is 



The Levy stable distribution has the following long-tailed behavior for large x 



~ a\x\ ^ X 



±oo. 



(6) 
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The Levy stable distribution reduces to a Gaussian for a = 2. Unlike the Gaussian distribution, 
the Levy stable distribution does not obey the usual central limit theorem. However, it obeys the 
Levy-Gnedenko generalized central limit theorem and is stable under addition 2^, 23]. 

An inverse power-law asymptotic behavior for the jump length distribution [Eq. ([7])], leads to 
a diverging mean square displacement (MSD) at all times: 



{x\t)) 



oo. 



(8) 



This is in contrast to Brownian motion where the MSD is finite and is proportional to time. 
This is because, the Brownian walker takes small steps at each interval of time. On the other 
hand, a Levy flier takes small steps interrupted by very long jumps (flights) in between. As a 
result, the variance in the step size is infinitely large 20(]. Levy flights are observed in a variety 
of phenomena. Search strategies of bacteria and various birds and animals, where an occasional 
long jump interspersed by short steps offers them significant advantages |28| . hopping on a 

polymer chain 12h14||. diffusion in micelles 15|], optical transport in a Levy glass [161 ]. and energy 
diffusion in single-molecule spectroscopy llTI] are a few examples. Levy flights have also been used 
to model encounters between different species (predator-prey) [2^ . 

Levy flights can be described by a generalization of the Fokker-Planck equation, referred to 
as the fractional Fokker-Planck equation (FFPE) The FFPE for Levy flight in the 

presence of an external potential V{x) is given by 



dP{x,t) 



-D 



dx"^ J 
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d V'{x) 
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In the above, D is the generalized diffusion constant with dimensions [D] = meter" sec~^, V{x) is 
the potential, 7 is the friction constant, and m is the mass of the particle. The solutions of the 
FFPE for free Levy flight and for Levy flight in the presence of a linear potential and a harmonic 
potential in the overdamped limit have been given by Jespersen et al. 
use the FFPE in the Fourier domain. For a free Levy flight, it is 



361 ]. In their paper, they 



dPfree{k,t) 

dt 
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which can be solved to get 

Pfreeik, t) = e-''^^^"' PfreeiK 0). (11) 

Imposing the initial condition Pfj.ee{x,0) = 6{x) leads to Pfreeif^^t) = e^^''^'"* which on Fourier 
transformation leads to the solution 

PM = ^ £ e-m:,-' = (^) . (12) 

For a Levy flight subjected to a linear potential of the form V{x) = —Fqx, the FFPE is 

^^'"^^''^-^ D\kr-^k^]pun{k,t) (13) 



dt \ 7m 

and the solution subject to the same initial condition as before is 



PliniKt) = e V 7m M 



-fo* .A _ 1 T ( ^ -ym 



This leads to 

PHn{^,t) = Pfree (x " = (^^",0 (^^^ J • (15) 

Finally, for Levy flight under a harmonic potential of the form V{x) = ^Xx"^, the FFPE is 

= -—k^P{k,t)-D\krP,,r{k,t). (16) 
at 7m ok 

One now has to employ the method of characteristics [s^ to obtain the solution 



As a result, 



A..(.,t) = exp(-^!^*^(l-e--/-)). (17) 



t) = "^_.x./.^x l ^-,0 f 77 ^ ^— — TT^l • (18) 



2?m7 (1 - e-°^*/7m) y ' \{Dm-f {l - e-'^^'/^--) /{aX)) 



The path integral approach has found extensive applications in the theory of Brownian motion. 
As Brownian motion is a special case of Levy flights considered here, one would expect path 
integrals to be very useful for Levy flights too. For example, a path integral approach for CTRW 
was discussed by Fredriech and Eule [37| and by Calvo et al. [3j| for fractional Brownian motion. 
However, though the path integral for a Levy processes has been given earlier, it has never been used 
for evaluating a probability distribution associated with these processes. An interesting extension 
of quantum mechanics was suggested by Laskin 3^ 4^ . It consists of modifying the Schroedinger 

(2 \ ~ 2 
— -^ ) in place of the usual operator He has used a 

path integral approach and has evaluated some of the associated path integrals. Our approach in 
the following is similar to that of Laskin, and can be used to evaluate Laskin-type path integrals. 
Essentially, we have the same kind of path integral, but our functionals are all real, as we are 
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concerned with real stochastic processes and not a generahzation of quantum mechanics. 

There has also been an enormous amount of interest in barrier crossing by a Brownian particle 



since the seminal work of Kramers (see the review by Hanggi et al. 4^). More recently, the escape 



of particles acted upon by fractional noise was investigated in Refs . l44l . l45l ]. Barrier crossing 



,l47|. 



problem for Levy flights has been discussed by a few authors [331 . |4£ 

In this paper we will use a Hamiltonian path integral approach for Levy flights. This type of 
path integrals were introduced in the context of quantum mechanics, long ago by Garrod [48|]. In 
these, the integration is over all paths in phase space, which are not continuous. If the action is 
quadratic in the momentum variable p, then the momentum path integral can be done easily to get 
an action containing only the position and its derivative, and one gets the more common integral 
over paths in position space. In the case of path integrals for Levy flights that we consider, the 
action expressed in phase space variables is not quadratic in p, but it has the term \p\"' , where 
a is a number (0 < a < 2) [See Eq. (j25p ]. Integrals over p then lead to symmetric Levy stable 
distribution, and hence one has a path integral involving a product of A^(— > oo) such distributions, 
which appears quite complex. It may be because of this that the path integrals for these processes 
has not been pursued in the literature. We show that an alternate approach in which, one first 
integrates over all the intermediate position variables and then integrates over the momentum 
variables, is quite powerful and allows us to evaluate several path integrals. 

Levy flights in the overdamped limit have been explored by other authors j2d. l36l. 1491] . Many of 
the results that we find in this limit, using path integrals, have been obtained by other methods. 
Our method, however, allows us to evaluate the propagator for Levy flight in a harmonic potential 
with a time-dependent force constant which has not been obtained by solving the FFPE. The 
underdamped limit for a free Levy flight has also been studied, but only to a limited extent by 
earlier authors. Some of the propagators in the underdamped limit for the special case of Levy noise 
with a = 1 (Cauchy noise) has been obtained by West and Seshadri 50(] and Garbaczewski and 
Olkiewicz 5]| . Srokowski [491 ] considers a case where the process is driven by an Ornstein-Uhlenbeck 
equivalent for the Levy noise and not by the Levy noise itself. Both additive and multiplicative 
noise are considered in his paper. However, in the underdamped limit, the propagator for only 
the free Levy flight is obtained, that too under certain conditions. Using our strategy, we are 
able to obtain the most general propagator in the underdamped limit for free Levy flight and for 
Levy flight in a linear and a harmonic potential. Our results reproduce the existing results under 
appropriate conditions. Further, the fact that one can evaluate path integral for a time-dependent 
harmonic potential exactly, has prompted us to study "semiclassical"-like approximations for path 
integrals, which can be very useful in the analysis of Kramers-like barrier crossing problem for Levy 
flights. We shall explore this in a forth coming publication 52l ]. 



II. THE HAMILTONIAN PATH INTEGRAL 

In this section we shall write a path integral expression for the propagator for Levy flight, 
subject to an arbitrary potential V{x). We start with the characteristic functional for a noise rj{t) 
defined by 

= (e*/o''<^*pW'?W). (19) 

The expectation value in Eq. (fT9|) may be written as a path integral over the noise variable, 'q{t), 
as 

V[p\ = [ Dr] /(T "^^f (20) 
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P[r]] is the probability density functional for rj{t). When the Levy noise is delta correlated, the 
most general form of the characteristic functional for Levy flights is [H^ 

= Jo dt {inp(t)-D\p(tT (l-i/3^^(p(t),a)) } ^21) 

where, a, /3, and uj{p{t),a) are as defined for Eq. ([3]). D plays the role of determining the 
noise intensity. Our notation needs a word of explanation. If the functional dependence is shown 
within usual brackets as in Pfree{x, t) in Eq. (jl2p . then it is just a function of the variables within 
the bracket. On the other hand, if we show the independent variable within square brackets as 
in VIp] of Eq. (j2ip . then this means that P is a functional of the function p{t). Also, note 
that P is used to denote the probability density functional, while V is used for its characteristic 
functional. We will describe our procedure for the simplest form of the characteristic functional with 
H = (3 = 0. However, our method will be applicable even for the general form of the characteristic 
functional. We will take D to have the dimensions [D] = meter°/sec, with [x] = meter and that of 
[p] = meter""^. P[r]] may be expressed in terms V[p] as 

= J Dp e-'^o<ii P(t)vit)^-Dj^fdt\pit)\'^_ (-22) 
We shall make use of the discretized version of the above integral, which may be written as 

— j / dpi dp2 / dpN exp f -At {ipnJln + P> IPnTI I (23) 

where we have divided the time interval (0, T) into A^(— > oo) equal intervals, each of length At 
and also, t„ = nAt and ry„ = f^" ^ dt r]{t). rjn ai^e independent Levy variables with the probability 

distribution for r?„ being given by P(r/„) = ^ dpn exp{-DAt - ipnr]nAt). 

The quantity that we are interested in is the probability density that a particle of mass m 
executing a Levy flight and starting at xq at the time t = would reach at the time t = T. The 
particle is assumed to obey the stochastic differential equation 

dx V'(x) , . 
dt rri'y 

rj(t) is assumed to be a white Levy process, having the characteristic functional Vlp] = 
e~^/o ^*lp(*)r^ which specifies all its properties. Using the above equation, and Eq. (j22p . we 
can write the probability density functional for x{t) as 

P[x] = I OpJe-^'i^"' \vi^Te-'^^''' (25) 

J is the Jacobian involved in the transformation from r]{t) to x{t). The propagator P{xf,T;xQ,0) 
can now be written as 

P(x^,r;xo,0) = r^^^^'' Dx [ DpJe-^Io^' W^T '^^ (26) 

Jx{0)=xo J 

At this point, it should be noted that the above equation resembles the Hamiltonian path integral 
of quantum mechanics, which was introduced long ago by Garrod [48|]. It is better if one also 
introduces the discretized version of the above path integral, so that it is possible to check for the 
correctness of the numerical factors, if necessary. It is convenient to take the discretized version of 
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Eq. (pij) to be given by 

^ (s^n— l) /■nrr\ 

X7 + = (27) 

At rwy 

With this choice of discretization, the Jacobian for the transformation J becomes unity, and the 
discretized version reads 

— / dxi.. dxN-1 / dpi.. dpN e-^^'^-=^\P"\" X 
^ *E„=iP.^^^ + ^^jAi^ ^^^^ 



Barkai a/. [5j] have derived a fractional generahzation of the Feynman-Kac equation for 
functionals of sub-diffusive continuous-time random walks. They have also derived a backward 
equation and a generalization to Levy flights. Solutions are presented for a wide number of appli- 
cations including the occupation time in half space and in an interval, the first passage time, the 
maximal displacement, and the hitting probability. For the particular class of processes that we 
are considering, it is quite easy to derive such an equation [Eq. ([9]) of this paper] using the above 
path integral. The details of this calculation is given in Appendix A. In the following, we shall 
write Eq. ^ as 



P(x/,r;xo,0) = / Dx Dp e-^Jo IfWI g-'Jo - yy^>y^^—)_ (29) 

ix(0)=xo J 

Note that we have adopted the Ito discretization, which is very convenient for the problem as 
there is no contribution to the Jacobian from the potential, V{x), in Eq. (|27p . The alternate 
Stratanovich version of discretization is 

H = Vn (30) 



At 2m7 



and results in a more complicated expression for the Jacobian [55 1. 

If a = 2, then the integrals over p„s in Eq. (|28|) are simple Gaussian integrals and can be done 
easily to get the usual path integrals for Brownian motion. For other values of a, each integral 
over pn will result in a stable distribution of the form 

'(31) 

Therefore, Eq. (|28p will become a product of N such stable distributions, where N ^ oo and 
this product has to be integrated over Xn, which seems formidable. This is probably the reason 
why a path integral approach for Levy flights has been used rarely in the literature. 

Our approach to the path integral in Eq. ([29|) is simple. If V{x) is at the most a quadratic 
function of x, then it is possible to first integrate over the position co-ordinates exactly leaving 
us with a product of Dirac delta functions involving the p„s. The integrals over p„ can then be 
performed easily to get the final answer. Using the procedure delineated here, we will obtain 
P{xf,T;xQ,0) for the free Levy flight and Levy flight in the presence of linear and harmonic 



potentials and thereby, reproduce the results of Jespersen et al. 36(]. For the harmonic potential. 



we will present results for both time- independent as well as time-dependent force constants. The 
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result for time-independent force constant is given in the paper by Jespersen et al. [36|]. However, 
the propagator for the Levy flight in a harmonic potential with a time-dependent force constant 
is a new result. Also, using the same approach, we will obtain the most general propagator in the 
underdamped limit for free Levy flight and for Levy flight in a linear potential and a harmonic 
potential with a time-independent force constant, which are also new results. The procedure can 
also be applied to find "semiclassical" approximations for more complicated V{x) and can be used 
to analyze the Kramers problem for a Levy particle SJ]. We adopt the following notations: In the 
overdamped limit, the propagator will have the superscript 'od' {P"'^), and in the underdapmed 
limit, it has the superscript 'ud' (P"*^). We also write p{T) = pT and p(0) = po- 



III. APPLICATION TO VARIOUS POTENTIALS IN THE OVERDAMPED LIMIT 

A. Free Levy flight 

For a free particle, V{x) = 0. For this case, Eq. (f28]l becomes 

(1 \^ /"°° /"°° JV 

— / dpi.. (ip^ e-^^*^«=ilP"l"e-'P^"/e*Pi"o X 

JJ / dxn e*(P"+i-P")^". (32) 



N-l 

n 

n=l " 

After performing integrals over all XjS, the resulting expression is 



P]t,{xf,T-x^,Q) = lim — / dpN e-^IP^l"^*e-^P'---^ TT / e"^l^"l'^^*<5(p„+i - p„) e^^^"». 

(33) 

Integrating over pn with n = 1 to — 1 gives 

Pfi,{xf,T;xo,0) = —j dpN e-^^l^'-l"e-*P-(^/-^«), (34) 



which leads to [36|] 



<e(x/,T;xo,0) = (^^) • (35) 

In the following, one can proceed with similar discretized integrals to evaluate the propagator 
for other problems. The discretized version of the path integral is tedious to perform, and it is 
far easier if one adopted a continuum version of the same integral. To illustrate the approach, we 
will do the above integral using the continuum approach. The approach has the difficulty that the 
propagator is determined only to within a multiplicative factor. However, this is not a problem 
as one can determine the factor using the normalization condition on the propagator. In all the 
further calculations, we will use the continuum version only. The propagator for free Levy flight 
in the continuum version is 



Pfi^{xf,T;xo,0)= Dx Dpe-~^^o'i^\p(^Te"'^o'itpii)^it). (36) 

Jx(0)=xo J 
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On performing the integral in the exponent, f^dt p{t)x{t), by parts, we get 

'■x{T)=Xf 



J J x{ 

= J Dp e^^^odt \pitTe-iiPTXf-poxo)^^p(^^^^_ (-37) 



Dx e* •'o 

x{0)=xo 



In the above equation, stands for Dirac delta functional, which results from the path integral 

over X. Dirac delta functional implies that 

p{t) = =^ p{t) = po, a constant. (38) 

On performing the integral over p{t), taking into account the delta functional and determining the 
multiplicative factor from the normalization condition, we obtain 

PUx„T;xo,0) = ^ lypo e--l-l-e--(^.-^o) = ^^L^, (^^) . (39) 
This is the same result that we obtained using the discretized version of the path integral in Eq. 



B. Linear potential 

We now analyze Levy flight under the influence of a linear potential of the form V{x) = —Fqx. 
The propagator is 



PlZ{xf,T; xo,0) = Dx Dp e'^Jo '^^ IPWI g-Jo ^v^v^.^;-;^; _ ^^q-^ 

Jx{0)=xo J 

On performing the integral in the exponent, Jq dt p{t) (x{t) — by parts. 



/T , \ Fn rT , ^ rxiT)=Xf „ 

Dp e-^/o"-^* \pitTe-i{PT-f-Poxo)^iJi /o dt p{t) / ^^^^dt pUMt) 

Jx{0)=xo 

(41) 

Integrating over all the paths in position space with fixed end points, we get 

Pitixf,T-xo,0) = j Dp e-^/o^'^t \p(t)re-i{PTXf-poxo)^i^!^dt (42) 

The Dirac delta functional (5[;j(t)], implies that 

p{t) = p{t) = Po, a constant. (43) 

Performing the path integration over p{t) accounting for the delta functional and determining the 
multiplicative factor such that P^"^ is normalized, gives 

Pit{xf,T-xo,0) = PfL {xf - FoT/mj,T;xo,0) . (44) 
If xq = 0, this is just the result of Jespersen et al. presented in Eq. ([T5]). 
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It may be noted that for /i 7^ and /3 = in Eq. ()2ip , the equations will turn out to be exactly 
the same as Levy flight in a linear potential. For the case where /? 7^ 0, this procedure, where we 
first integrate over position coordinates first, is still valid. As an example, let us consider free Levy 
flight where /? 7^ and 7^ 0. After performing the integration over position variables, we will be 
left with the integral 



1 



00 



The final distribution is Levy stable, just that it is not symmetric, as expected. 



C. Harmonic potential with a time-dependent Force Constant 



We will now consider Levy flight in a harmonic potential with a time-dependent force constant, 
for which the results have not been reported in the literature. It is interesting to note that our 
method is very easy when applied to this problem, while in comparison solving the corresponding 
fractional differential equation would be much more involved. We write the potential as V{x) = 
^A(t)x^ and find an analytical expression for the propagator using the Hamiltonian path integral 
approach. The path integral for the propagator is 



P^Kx/,T; xo,0) 



x{T)=Xf 



Dx 



Dp e 



lx{0)=XQ 

On performing the integral in the exponent, Jq dt p(t)x{t), by parts 

rxiT)=Xf 
'x{0)=xo 

Dp ^0 "^^ lp(*)l"g-«(PTa;/-poa::o)j 



(46) 



! 



m - ^p{t) 



m7 



The delta functional leads to the relation 

m 



Pit) 



m7 



p{t) = ^ p{t) = poe^o dt' A(t')M7. 



(47) 



(48) 



After performing the integration over the delta functional and making a change of integration 
variable pq = po e^o -^(*)/n^7^ ^j^e propagator, P^^^(x/, T; xq, 0) along with the multiplicative 
factor is 



P;:i{xf,T;xo) 



[Dj;^dt e-/.^^*' ^i^/^^y^'" 1^ [Dj;[dt e-^f^^'' ^(*')/-^) 



. fT at' ^ 

■J m7 



l/a 



(49) 



For time-independent force constant, the result in Eq. (|49p reduces to the following 

/ 



P^l,(xy,T;xo) 



1 / -^a.yj 

(Z)m7(l — e ""T )/(aA)j W Z)m7(l — e )/(aA)j 



AT 

Xf — XqC 



l/a 



(50) 
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This result is identical to that of Jespersen et al. in Eq. (jlSp for the initial condition xq = 0. Clearly 
as T — )• oo the above function approaches /„ o ( -, — -tt777 ) > showing that the probability 

{Dm-yj ' ' \ (Dm-y) ' ) 

distribution attains a steady value at infinite time. The average energy of the particle is given by 
the expression |A(a;j). Since the mean square displacement of a Levy distribution always diverges 
[Eq. ([8])], the average energy of the particle also diverges. 

IV. LEVY FLIGHT IN THE UNDERDAMPED LIMIT 

The use of our procedure, enables us to obtain the most general propagator in the underdamped 
limit, which is a new result. In the underdamped limit, the inertial term in the Langevin equation 
cannot be ignored and the motion is governed by the equation 

x(t) . , , V (x(t)) , , 

-^+x{t) + — ^-±^=r]{t). 51 

7 m7 

Here also, we will use the Hamiltonian path integral formulation in its continuum version. We will 
present results for the free Levy flight and also in the presence of a linear potential and a harmonic 
potential with a time- independent force constant. 



A. Free Levy flight 



The propagator for free Levy flight in the underdamped regime is 

rxiT)=Xf [■ ^ _ ^ ^^^^^ _. ^T^, „u\(m. 

'x(0)=xo 



PjUxf,T-xo.^) = r^^^^'' Dx [ Dp e-^fo'^ mi^-^Io'^^ P(t)(^+^W)_ (53) 

Jx(0)=xn J 



Integrating dt p{t) {^^^ + x{t)^ in the exponent by parts we get 

V 7 / 7 

\t (^P^-pit)^x{t), (53) 

where, vq and Vf are the initial and final velocities, respectively. Since the differential equation 
governing the motion is second order [Eq. (j5ip ]. the propagator will depend on the initial and final 
positions as well as the velocities. The propagator is 

Dp e-^fo^' +PT^f-Po^o) ^ 

Dx e"'^^"' 

Dp e-^/o \pitTeA— ^ [p(t) - ^pit)] . 

Dirac delta functional leads to the condition 

Pit) - im = 0, (54) 
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which has the solution 



.\p{t) = ci + C2e 



7* 



(55) 



where ci = ^^^~^°t and C2 = ^^~^t • On performing the integral over the delta functional, the 
path integral over all ps reduces to just two integrals over px and po- Therefore, the unnormalized 
propagator is 



Pfree{Xf,Vf,T;Xo,Vo,0) 



dpQ 



dpT e 



exp 



I < PT 



Vf vq poe'^'^ - PT 

— -P0— + — Zt 

7 7 e ' 



PT f s 
— {Xf-Xo) 



(56) 



After normalizing, the propagator becomes 



exp 



dpo / dpT e 



-Dfodt 



po(i-'=~^')+pt('=~""-'=~^'^) 



f , Po - Pre 
Mpt Po— + — j ZZt 

7 7 1 — e 



7T 



-(x/ - Xo) 



(57) 



For a = 2, the above expression corresponds to Brownian motion, the integrals can be performed 
exactly and, the propagator we obtain matches with the result given in the seminal paper by 
Chandrasekhar 56(]. For a 7^ 2, it seems difficult to evaluate the double integral in Eq. ()57p . In 
order to simplify the evaluation of the propagator, we make a change of variables from {po,Pt} to 
{gi, (72}, such that px = Qi and po = qiq2- The details of this transformation are given in Appendix 
B. On making the change of variables, we get 



Pfreeixf,Vf,T;xo,vo,0) = ^^^^^ _ ^^^^^ J dq2 dqi gie"®il''il°cos(Aigi) 



(58) 



where, 



D 



T 



dt 



e-7i _ e-7T 



-It 



1 



-7T 



+ 92 



1 



-7T 



Ai = ^ 



V e ^'^{xf - xq) 



Q2 — 



[Xf - Xo, 



7 



1 



7 1 - e-^^ 

The integral over qi can be performed exactly for a = 1 and written as 

1 



7T 



ryud 

free 



'xf,Vf,T;xo,vo,0) 



27r27 (1 - e-T^) 



02 ^2 

dq2 o- 

-00 {Qj + Ajf 



(59) 



(60) 



For other values of a, the integral over qi cannot be performed exactly. However, we expand the 
cosine function in Eq. (j58p and then perform the integral over qi, which results in 



PfLixf,Vf,T; xo,vo) = ^ 



1 



7Q (1 - e-T"^) 



n=0 



(2n)! 



a 



Af 



®1 



2(n+l)/a ■ 



(61) 



The sum does not converge for a < 1. As T — >■ cx), Pj^^^{x f,Vf,T; xo,vo) ^ because 0i — >• 00 in 
this limit for any non-zero value of q2- 
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More specific propagators can be derived from Pj^^^{x j , v f ,T; xq, vo,0) . For the case of Brow- 
nian motion, the propagators Pj^^^{x f,T; xo,Vo,0) and PJ^^^{vf,T;vo,0) were given in Chan- 
drasekhar's paper. Similar exact results can be obtained for Levy flights using our formalism. 
The details of the calculation are provided in Appendix C. The propagator 

/oo 
dvf Pp,,{xf,Vf,T;xo,vo,0) 
-oo 

La,o (^] , (62) 



'2 



with 



and 



fl-e-^'Y I? B (1 - e-T^; 1 + a, 0) , , 



In the limit T oo, ©2 oo and therefore, in this limit, Pjreei^f^ ^0; "^0; 0) — ?• 0. 

fr 



On integrating over all possible values of of PV'^^g{v f , x f ,T; vo,0) , we get 



/oo 
dxf Pp,,{xf,Vf,T;xo,vo,0) 
-oo 



where 



and 



03 = ^ (1 - , (66) 



= Vf — vqc "''^ . (67) 
The stationary distribution for PJ!^^^{vf,T;xQ,VQ,^) is obtained by letting T — >• oo, which gives 



This reduces to the Maxwell Boltzmann velocity distribution for a = 2. PJ^^^{xf,T;xQ,VQ,0) and 
Pi^eJ:!Jf,T;vo,0) for the special case of Cauchy noise (i.e. a = 1) have been obtained earlier 
[sol.lsil. For Q = 2, we reproduce the expressions given in Chandrasekhar's paper [sH ]. 
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B. Linear Potential 



The propagator for the Levy flight under a linear potential of the form V{x) = —Fqx in the 
under damped regime is 



Pli^{xf,Vf,T;xo,vo,0)= / Dx / Dp e- 



7771 



(69) 



Since, V'{x) is independent of the position in this case, the part of the above expression dependent 
on position is exactly similar to the case of free Levy flight. On following the same procedure as 
we did for free Levy flight, we obtain the normalized propagator to be 



P^i{xf,Vf,T-XQ,v^,{)) 



1 



47r27 (1 - e-T^) 



00 roo 

dpo / dpT e 

-00 -/ —CO 



exp 



-I { PT - PO + 



7 



7 



1 - e- 



FqT 

m7 



(70) 



From the above expression it is easy to see that the case of Levy flight under a linear potential is 
equivalent to free Levy flight with Vf -^^ {vf - , vq (^vq - ^) and xj ^ (^xj - . On 
making these substitutions in the propagators for free Levy flight we can obtain the corresponding 
propagators for Levy flight in a linear potential. 



C. Harmonic Potential 



In the underdamped limit of friction, we analyze Levy flight in a harmonic potential of the form 
V{x) = iAx^, with a time-independent force constant. The propagator, P^^^{xf,Vf,T;xo,vo), is 



Pl^S"ixf,Vf,T;xo,vo)= Dx Dpe 



Xx{t) 
7717 



Integrating dt p(t){^^ + x{t) + ^^^} in the exponent by parts gives 

\tpit)(^+±it) + ^] = PTVf-PO-O-PTXf+PQXo 



7 



rwy J 



7 



+ PTXf -POXO + 



^ dt (^-p{t) + —p{t)]x{t). 







7 



Using this result in the propagator we obtain, 
Pharixf,Vf,T;xo,vo) = J Dp e 



rj, .f PT^f-PO^O-PT^f+P'o^O , ^ 
D ff^dt \p{t)\"^\ J +PTXf-poXOj 



Dp e-^^o'it 



+PTXf-poXo 1 



Pit) - ^p{t) + -p{t) 

m 



(71) 



(72) 



(73) 



(74) 
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Delta functional implies, 



Pit) - im + -p{t) = 



m 



■'■ P{t) = e 2 ci cosh [ —t ] + C2 sinh —t 



(75) 
(76) 



where /3 = -1/7^ — 4—, ci = po and C2 = ^-^-^^ — ^/"fl^t^^ ^^'^ ■ Though we have shown the working 



only for the case where 7^ > ^) our analysis is equally valid in general. Integral over the delta 
functional will result in the in a double integral over px and pq. The resultant unnormalized 
propagator is 



Phari^f^V.fiT;Xo,Vo) 



-Dj^dt 



dpo / dpT e 



7t PT <! sinh ( ) +P0 si^h ( '^^'^ *^ 



PT^-PO^ + ^ ^ % ^ '-Xf ^ ^ ^ 



(77) 



Upon determining the normalization constant, the propagator is 



Phar{xf,Vf,T]Xo,Vo) 



7T 
e— /3 



dpo / dpT e 



p„ / 7-/3coth(^) _7T 



/5-(^)^)+^(^-.0-^i^^^!^xo+e^/5c.K(^)^) 



In order to obtain a more elegant expression for the propagator, we perform a simple change of 

7T 



variables where pt = — and po 



7 sinh 



and obtain 



^ rco roc 

Pharixf,Vf,T;xo,vo) = ^ J dpQ J dpT e 



_7t 



■ye T pQ sinh(;|^)+pT 



7-/3cothf 7T 



x^+e-T-^ csch(^)^ |+po( -voe-^ smh(^) + 



Pxf 7sinh(4^)+/3cosh(4^) _ fT 



e i xo 



It is possible to get the stationary distribution by allowing T —?■ 00. Noting that we are considering 
the case 7 > /3, and letting T ^ 00 leads to 



Phar,stixf,Vf) = ^ j j dPT e 

-«{pt(^'/ + ^^7)+P0^} 



7e-^(posmh(4^)+pTe-'"/2) 



(80) 



The above equation implies that P^ar st^'^'^f^'^f) two-dimensional Fourier transform of a well 

behaved function. This guarantees that the steady state distribution is well-behaved for any 
2 ^ Q > 0. 
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In Eq. (j79p . after performing the coordinate transformation from {po,pt} to {(?i,(?2} as we did 
for the free particle, we get 



where, 



47r2 



a 



-oo Jo 



COS ixiqi) 



T 

ni = D I dt 
'o 



76 



sinh ( ^^"^ 
-7t/2 I V 2 



/™x + 92 sinh 
sinh I ^ ^ 



Xi = I ^'z + ^ a;/ + e 2/3 csch ( ^ 1 y 



7T 



+<?2 —"Woe 2 sinh 



(cosh ( ^ 



V 2 y 



_2L 

e 2 xo 



(81) 



(82) 



For Q = 1, the integral over qi can be performed exactly resulting in the expression 



/3 



Pht {xf,Vf,T;xo,vo) = / (ii?2 



oo (Of + xl) 



(83) 



This enables us to calculate the time evolution of the probability distribution easily (see Fig. 1). 
For other values of a, for which the integral over qi cannot be performed exactly, we expand the 
cosine function in Eq. (|8ip as a sum and perform the integral over qi. This results in 



PhariXf,yf,T;Xo,Vo) 



(3 



47r2 



a 



n=0 



(-1)" ^/ 2(n+l) 
(2n)! V " 



X?" 



2{n+l)/a ■ 



n 



(84) 



This sum does not converge for a < 1. 

Specific forms of the propagator PJ^^^{xf,T;xo,vo,0) and P^^^(f/, T; xq, fo, 0) can be obtained 
from P^^^{xf,Vf,T]XQ,VQ,Q) by integrating over Vf and Xf respectively. P^^^(xj, T; xq, fo, 0) is 



P]ii{xj,T-XQ,VQ,Q) 



dvf Pha,^{xf,Vf,T;xo,vo,0) 



1 r / ^2 



1/a' 



a 



l/a 



(85) 



where, 



n2 = D I dt i y e-^*/2 sinh f 



/5. 



V 2 V/3 



1. ff3'P\ 7 . , f(3T\\ ^ 2vo ^ . ^ (PT\ 
X2 = Xf-[ cosh I — I + ^ smn ( — ) I ^ ^ xq - — e 2 smh ( — I . 

The stationary distribution for the propagator P^^^(x/, T; xq, 0) is 

dvf Pr,i{xj,Vf,T-x,,v,,Q) = ^L^,o[^ji] 



Phar,st{Xf) 



(86) 



(87) 
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where, 



X2,st = Xf. 



(88) 



Also, 



P;:^,{vf,T;xo,vo,0) 



dXf Pharixf,Vf,T;Xo,Vo,0) 
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1/a 



L. 



a,0 



X3 



a 



l/a 



(89) 



where. 



= D dt 



T 



-ye .'/-(cosh(^*)-^smh(^|t 
,3 . , I j.„, (f) - eosh (f )) e-^^/'^.„ . ( '^' - ^^'J'" I e-^/^.„.(.0) 



The stationary distribution for the propagator P^^^(x/, T; xq, "Wo, 0) is 



har,st 



/CO 
dXf Phar {Xf ,Vf,T; Xo,Vo,0) 
-oo 



2n 



P J. I X3,st 



l/a' 
3,st 



l/a 
3,st / 



(91) 



where. 



^3,st = D dt 
Jo 

X3,st = Vf. 



76 7^2 (cosh -^sinh^^t 



(92) 



For Q = 2, we obtain the results of Chandrasekhar [5f 

The propagators in the case where 7^ < ^ can be obtained from the expressions derived for 

7^ > ^ by replacing cosh[^] with cos[^] and sinh[^] with isin[^] where, (3 = \J ^ ~ 1^ ■ 
For 7^ = ^ case, the appropriate propagators can be obtained by replacing cosh[^] with 1 and 
sinh[^] with ^ in the expressions obtained for 7^ > The expectation value of the energy of 
a particle executing Levy flight in a harmonic potential in the underdamped limit will be given by 
^'m{v'j) + ^A(xj). Since both the position and the velocity distributions are Levy stable, as given 
in Eq. (j85|) and ([89]) respectively, both (x^) and (v'j) diverge. 

Phase space distributions for Levy flights under a harmonic potential have been explored previ- 
ously (see references 57, 5^). In the case usual Brownian motion, Xf and Vf are uncorrelated and 
one obtains elliptical contours for the phase space distribution. However, for Levy flights this is no 
longer true - an inhomogeneous stationary phase space distribution was observed 57, [HI]. Since, 
our method gives exact expressions for propagators at any time, we can study the time evolution 
of the phase space distribution, while the previous papers were concerned only with the steady 
state distribution function. In Fig. (1(a), 1(b), 1(c), 1(d)), we show the time evolution of the phase 
space distribution for a = 1 with the initial conditions, Vf = and Xf = 0. In Fig. (1(e)), we 



also show the stationary distribution calculated simply as the product of P}^^j,{xf,T; vq,xq,0) and 



PJ^^^{vf,T;v(),XQ,0) in the large T limit. We observe that there is correlation between xj and Vf 
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which does not die off as T — )• oo. Our fin ding s are consistent with the inhomogeneity in phase 



space distribution reported in the literature [57|, [5 



V. CONCLUSIONS 

We have shown that Hamiltonian path integrals offer an appealing way to study Levy flights 
in both the overdamped and underdamped limits of friction. Though some results are already 
available, the method adopted to obtain them is interesting. In addition, they enable us to obtain 
expressions for several propagators for Levy flights, which seem difficult to obtain using other 
methods. 
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VII. APPENDIX 
A. Fractional Fokker-Planck Equation from the path integral 

Here we derive the FFPE starting from the discretized of the path integral as given in Eq. ()28p . 
Let P{x, t + Ai; xq, 0) be the probability of finding the particle at x at the time t + lS.t given that 
it started at xq at the time t = 0. Then the path integral prescription of Eq. (j28p enables us to 
calculate this probability density in terms of P{x' , t; xq, 0), the probability distribution at the time 
t, as 

1 roo roc ( v'(x') . ^ 

P{x,t + At;xo,0) = / dp dx' e-^l^'l"^Vn(^-^)+— Mp(:E',t;xo,0). (93) 



2tt 



Since At is infinitesimally small, we can expand the exponential ev t"" / up to linear 

order in At and obtain, 

1 roo roo / V'(r')\ 

P(x, t + At; xo, 0) = — / dp dx' (l- DhrAt + ipAt^^ e^^'^^"^ ^P(x', t; xo,0). 

(94) 

Defining the operator I — j by 

(o2 \ <^/2 1 roo roo 

--^ ^0' 0) = ^ / '^P dx'D\pre'P<^^-^'^P{x', t; xq, 0), (95) 



we can write Eq. (|94p as 

Pfrr f- -\- At \ -D I - 

dx"^ J dx 7m 



q2 \"/2 Q yi(x\\ 

P(x,t + At;xo,0) = P(x,t;xo,0) + Atl +- lp(x,t;xo,0). (96) 
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From the above equation, we obtain the FFPE given in Eq. ([9]) which is 
aP{x,t;xo,0) ^ , B V\x) 



B. Transformation of coordinates from {potPt} to {qi,(?2} 

The propagators PJ^^^{x f , v f ,T; xq, vq,0) and PJ^I^^{xf,Vf,T;xo,vo,0) can be written in a gen- 
eralized form as 

-OO J — oo 

which shows the formal dependence of po and pT- We will not work out the steps explicitly 
for Pj^^g {xf,Vf,T;xo,vo,0) and P^f^ (x/ , t^/ , T; xq , I'd , 0) but only for P^en • The propagators of 
interest can be obtained from Pgen- 

Using the transformation {po,Pt} to {qi, 92} with pT = qi and 

Po = Q1Q2, we get 

/oo /-oo J, 

dq2 / dgi -Zie-l^il^/o Ht)+mi.\- ^os{{g + hq2)qi). (99) 
-00 ^0 

Upon expanding the cosine series and performing the integral over gi, we get 



n=0 



/g [a(f) + 6(t)g2h 



Note that the above sum can be evaluated exactly for a = 1 or 2 and lead to Cauchy or exponential 
functions. Further, the sum does not converge if a < 1. 

C. Calculation of the propagator PJ^^^{xf,T;xo,vo,0) 

Our procedure enables us to calculate the most general propagator in the underdamped 
regime for both free Levy flight and that under a linear and a harmonic potential. From these, 
we showed results for more specific propagators such as PJ^^^{xf,T;xo,vo,0), PJ'^^^{vf,T;vQ,0), 
PJ^^^{xf,T;xo,VQ,0), and Pl^^j.{vf,T;xQ,VQ,0). Here, in the appendix, we will show the steps in- 
volved in calculating Pf^ee^'^f^'^'^ xq,vq,{)). The rest of the propagators can be obtained in a similar 
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fashion. 



/oo 
dVf PfreeiXf,Vf,T; Xo,Vo, 0) 
-oo 



oo 



AiT^-f (1 - e-T^) 



dpo / dpT e 



3(l-e-Tt)-p^(e-7T„,-7t) 



CO 
7T 



dpQ / dpT e 



27r(l-e-^^) 



l-e-7'J' 



(101) 



Performing the integration over the delta function we obtain 



1 



27r (1 - e-^^) 



dpQ e 



(1-e- 



a,0 



1/a 



.(102) 



[1] H. Scher and M. Lax, Phys. Rev. B 7, 4491 (1973). 

[2] Q. Gu, A. Schiff, S. Grebner, and R. Schwartz, Phys. Rev. Lett. 76, 3196 (1996). 

[3] H.-P. Muller, R. Kimmich, and J. Weis, Phys. Rev. E 54, 5278 (1996). 

[4] A. Klemm, H.-P. Miiller, and R. Kimmich, Pliys. Rev. E 55, 4413 (1997). 

[5] F. Amblard, A. C. Maggs, B. Yurke, A. N. PargcUis, and S. Lciblcr, Phys. Rev. Lett. 77, 4470 (1996). 

[6] G. M. Viswanathan, V. Afanasyev, S. V. Buldyrev, E. J. Murphy, P. A. Prince, and H. E. Stanley, 

Nature 381, 413 (1996). 

[7] A. M. Reynolds and M. A. Frye, PLoS ONE 2(4), e354 (2007). 

[8] R. P. D. Atkinson, C. J. Rhodes, D. W. Macdonald, and R. M. Anderson, OIKOS 98, 134 (2002). 

[9] G. Ramos-Fernandez, J. L. Mateos, O. Miramontes, G. Cocho, H. Larralde, and B. A.-O. B, Behavioral 
Ecology and Sociobiology 55, 223 (2004). 

[10] F. Bartemius, Fractals 15, 151 (2007). 

[11] F. Bartumeus, J. Catalan, U. L. Fulco, M. L. Lyra, and G. M. Viswanathan, Phys. Rev. Lett. 88, 
097901 (2002). 

[12] 1. M. Sokolov, J. Mai, and A. Blumen, Phys. Rev. Lett. 79, 857 (1997). 

[13] D. Brockmann and T. Geiscl, Phys. Rev. Lett. 91, 048303 (2003). 

[14] M. A. Lomhoh, T. Ambjornsson, and R. Metzler, Phys. Rev. Lett. 95, 260603 (2005). 

[15] J. Klaftcr, A. Blumen, G. Zumofen, and M. F. Shlesinger, Physica A 168, 637 (1990). 

[16] P. Barthclcmy, J. Bertolottil, and D. S. Wiersma, Nature 453, 495 (2008). 

[17] G. Zumofen and J. Klafter, Chem. Phys. Lett. 219, 303 (1994). 

[18] E. W. MontroU and G. H. Weiss, J. Math. Phys. 10, 753 (1969). 

[19] E. W. MontroU and H. J. Scher, J. Stat. Phys. 9, 101 (1973). 

[20] R. Metzler and J. Klafter, Physics Reports 339, 1 (2000). 

[21] J. Klafter, A. Blumen, and M. F. Schlesinger, Phys. Rev. A 35, 3081 (1987). 

[22] P. Levy, Calcul des probabilitcs (Gauthier-Villars, Paris, 1925). 

[23] G. M. Zaslavsky and U. Frisch, in Lecture notes in physics. Vol. 450, edited by M. F. Schlesinger 
(Springer Vcrlag, Berlin, 1994). 



20 



[24] A. V. Chechkin, R. Metzler, Klafter, and V. Y. Gonchar, in 

Anomalous Transport: Foundations and Applications , edited by G. R. R. Klages and I. M. Sokolov 

(Wiley- VCH Vcrlag, Wcinhcim, Germany, 2008). 
[25] R. Metzler, A. Chechkin, and J. Klafter, Encyclopaedia of Complexity and System 

Science, Levy statistics and anomalous transport: Levy flights and subdiffusion, (Springer- 

Science+BuisinessMedia, LLC, 2009). 
[26] B. V. Gncdenko and A. N. Kolmogorov, Limit Distributions for Sums of Independent Random Variables 

(Adison-Weasley, Reading, 1954). 
[27] W. Feller, An Introduction to Probability Theory and its Applications, vols I and II (Wiley, New 

York, 1971). 

[28] G. M. Viswanathan, M. G. E. da Luz, E. P. Raposo, and H. E. Stanley, 
The Physics of Foraging: An Introduction to Random Searches and Biological Encounters (Cam- 
bridge University Press, 2011). 

[29] I. Goncharenko and A. Gopinathan, Physical Review Letters 105 (2010), URL 
jhttp : //dx . doi . org/10 . 1103/PhysRevLett . 105 . 190601 

[30] R. Metzler, E. Barkai, and J. Klafter, Europhys. Lett. 46, 431 (1999). 

[31] R. Metzler, E. Barkai, and J. Klafter, Phys. Rev. E 61, 132 (2000). 

[32] R. Metzler, E. Barkai, and J. Klafter, Phys. Rev. Lett. 82, 3563 (1999). 

[33] A. Dubkov, A. La Cognata, and B. Spagnolo, Journal of Statistical Mechanics: Theory and Experiment 
P01002, 1 (2009). 

[34] I. Calvo, R. Sanchez, and B. A. Carreras, J. Phys. A: Math. Theor. 42, 055003 (2009). 

[35] W. Ebeling and I. M. Sokolov, Statistical Thermodynamics and Stochastic Theory of Nonlinear Systems Far from Equil: 

(World Scientific Publishing Co. Pte. Ltd., 2005). 
[36] S. Jespersen, R. Metzler, and H. C. Fogedby, Phys. Rev. E 59, 2736 (1999). 
[37] R. Friedrich and S. Eule, arXiv:1110.5771vl (2011). 
[38] N. Laskin, Physical Review. E 62, 3135 (2000). 

[39] N. Laskin, in Fractional Dynamics: Recent Advances , edited by J. Klafter, S. C. Lim, and R. Metzler 

(World Scientific, 2012), pp. 393-428. 
[40] N. Laskin, Phys. Rev. E 66, 056108 (2002). 
[41] N. Laskin, Phys. Lett. A 268, 298 (2000). 
[42] N. Laskin, Chaos 10, 780 (2000). 

[43] N. Laskin, Commun. Nonlinear Sci. Numer. Simul. 12, 2 (2007). 

[44] P. Hanggi, P. Talkner, and M. Borkovec, Rev. Mod. Phys. 62, 251 (1990). 

[45] S. Chaudhury, D. Chattcrjec, and B. Cherayil, The Journal of Chemical Physics 129 (2008), URL 

pittp: //dx . doi . org/10 . 1 063/1.2969767. 
[46] A. Chechkin, O. Shusarcnko, R. Metzler, and J. Klafter, Physical Review E 75 (2007), URL 

|http : //dx ■ doi . org/10 . 1 103/PhysRevE . 75 . 041101( 
[47] B. Dybiec, E. Gudowska-Nowak, and P. Hanggi, Physical Review E 75 (2007), URL 

|http : //dx ■ doi ■ org/10 . 1 103/Phys RevE . 75 . 021109 
[48] CrCarrod, Rev. Mod. Phys. 38, 48371966). 
[49] T. Srokowski, Phys. Rev. E 85, 021118 (2012). 
[50] B. J. West and V. Seshadri, Physica 113A, 203 (1982). 
[51] P. Garbaczewskia and R. Olkiewicz, J. Math. Phys. 41, 6843 (2000). 
[52] D. Janakiraman and K. L. Sebastian, to be published. 

[53] A. Chechkin, R. Metzler, J. Klafter, and V. Gonchar, Anomalous Transport pp. 129-162 (2008). 
[54] S. Carmi, L. Turgeman, and E. Barkai, J Stat. Phys. 141, 1071 (2010). 
[55] A. J. McKane, H. C. Luckock, and A. J. Bray, Phys. Rev. A 41, 644 (1990). 
[56] S. Chandrasekhar, Rev. Mod. Phys. 15, 1 (1943). 

[57] I. M. Sokolov, W. Ebeling, and B. Dybiec, Phys. Rev. E 83, 041118 (2011). 
[58] Z. Cao, Y.-F. Wang, and H.-G. Luo, Phys. Rev. E 85, 042101 (2012). 



21 



(a) P^t{vf,Xf, T; 0, 0, 0) at T=l (b) P^"<f.(«/, xj,T- 0, 0, 0) at T=3 
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(c) P;^^,ivf,Xf,T;0,0,0) at T=5 (d) P^^,{vf,Xf,T;0,0,0) at long time 
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(e) Product of the Stationary 
distributions P^ti^f) and P^^tivf) 
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FIG. 1: The figures (a)-(d) show the time evolution of PJ^I^^{v j , x f ,T; vq, xo,0) starting with 
f = and xq = 0. Figure (e) shows the long time behavior of the product of P^^r^x f ,T] vq,xq, 0) 
and Pl^i{vf,T;v^,xo,Q). (a) P]ii{vf,Xf,T;Q,Q,Q) at T=l; (b) P^i{vf,Xf,T-QM at T=3; (c) 

-fw(^/'^/'^;0'0'0) at T=5; (d) (a) P^^^.{vf,Xf,T- at long time; (e) Product of the 
Stationary distributions Phari^f) ^^'^ ^haA'^f)- 



